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ABSTRACT 



A military aircraft deployment problem from the United States Transportation 
Command is modeled as a generalized transportation problem with side constraints and 
solved using a general purpose linear programming package. The model involves the 
assignment of military units and material to aircraft and the assignment of aircraft to 
missions in order to appraise the utility and to determine the assets required for prelim- 
inary military operation plans. A transformation of this model which aggregates vari- 
ables relating to the early or late delivery of requirements is also described. A specialized 
algorithm which separates an instance of the model into subgroups of independent time 
windows, finding the globally optimal solution by solving independent subproblems is 
also explored. Lastly, an integer rounding model is described which converts continuous 
solutions to integer in order to facilitate implementation of the former models with an 
existing post-solution processor. 

Excellent quality solutions are provided for problems involving nine routes, SO 
movement requirements distributed across two cargo classes invohing 200.000 short 
tons of freight, and 250 aircraft using four different aircraft types for each of 12 time 
periods. The problem, which has the potential of having over 10,000 variables, is re- 
duced significantly using variable reduction and the aggregation transformation. The 
reduced problem requires approximately 1.000 variables and 300 constraints and sol- 
utions are obtainable in under 14 minutes using the General Algebraic Modeling System 
on an S02S6-based personal computer. 
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I. INTRODUCTION 



In the modern world of rapid mobility, the difference between a successful and an 
unsuccessful deployment of military units and equipment will depend on the efficient 
allocation of transportation assets to the requirements which need to be moved within 
prescribed time intervals. These transportation assets consist primarily of military 
transport aircraft, ocean-going vessels, and rail and truck vehicles. The movement re- 
quirements may consist of any material item or person requiring transportation from one 
location to another such as an infantry division, petroleum products, heavy equipment, 
artillery, subsistence items, and many others. Small differences in the timing of deliver.’ 
or the effective use of assets could determine the difference between an action being won 
or lost [Ref. 1]. 

This thesis considers the question of whether a potential military operation can be 
initiated and sustained from the perspective of air transportation. It is valuable to know 
what assets will be used, where they will be needed and whether there are bottlenecks in 
the transportation process which significantly affect the operation. It is not necessary 
to develop detailed schedules and load plans for aircraft and movement requirements for 
this initial planning: creating such schedules and plans is beyond the scope of this work. 
However, the initial planning process for the allocation of transportation assets will need 
to be sensitive to the timing of delivery of these movement requirements and to the ef- 
fective use of the transportation assets. The problem of determining the appropriate 
deployment of aircraft given a set of requirements which must be moved within pre- 
scribed time intervals and with specific delivery times which are the most opportune for 
the objectives of the operation will henceforth be called the deployment problem. Past 
model development has shown that the complexities involved with combining different 
modes of transportation in a single model require solving an extremely large problem. 
Consequently, a modular approach which separates the air, land, and sea allocation 
problems is usually employed [Ref. 2]. 

The United States Transportation Command (USTRANSCOM), is the Department 
of Defense's agency which has been most concerned with the deployment problem. This 
organization was created to consolidate the global air, land, and sea transportation as- 
sets and requirements which are used by the military to meet national security objectives. 
While tasked with coordinating the development of deployment plans involving 
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transportation assets, one of USTRANSCOM's principal endeavors is to provide 
automated data processing systems to streamline the planning process. The Joint 
Deployment Agency (JDA). the predecessor organization of USTRANSCOM, 
recognized the potential of automated data processing to assist with mobilization plan- 
ning. The Joint Deployment System (JDS), a comprehensive computer planning system 
initiated under the auspices of the JDA, was designed to plan, execute, and monitor the 
deployment of units during both peacetime and times of crisis. Many reseach initiatives 
directly related to the deployment problem have been undertaken to support this system. 
[Ref 3] 

The first linear programming based effort at modeling and solving the deployment 
problem was the System for Closure Optimization Planning and Evaluation (SCOPE), 
a submodule of the JDS [Ref. 4], This system used a Benders decomposition method to 
solve a continuous variable model consisting of both sea and air assets over discrete time 
periods. The formulation separated the model into a master problem involving the as- 
signment of assets to delivery missions and a subproblcm involving the assignment of 
movement requirements to transportation assets. Because the solution method con- 
verged slowly or not at all. and because it did not model sea assets integrally, the Benders 
decomposition approach was considered to be unsuitable for the JDS system [Ref. 5 : 

p. 20]. 

The deployment problem for the allocation of sea assets has been studied by many 
researchers since the original SCOPE model. A continuous variable deployment model 
for both sea and air assets similar to that of the SCOPE model was undertaken by 
Collier [Ref. 5] using linear programming and a variable reduction heuristic. This model 
was developed to assist planners with designing feasible deployment plans and selecting 
appropriate modes of transportation for movement requirements. An alternative model 
for the deployment problem involving only sea transportation assets was developed by 
Lally [Ref. 6] using integer programming. This model avoided the integrality problems 
which occurred for sea assets in the SCOPE model. While both of these studies dem- 
onstrated the ability to use optimization techniques for moderately sized instances of 
problems, neither of these techniques have been employed in the JDS system. In its 
current form, the JDS uses simulation-based heuristics which are capable of providing 
solutions to highly detailed and very large resource allocation models. However, the 
JDS still has difficulty in providing solutions to large instances of problems in a timely 
manner [Ref. 7], and the solutions obtained cannot be considered to be optimal with 
respect to any measure of effectiveness. More promising developments involving opti- 
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mization techniques for solving large-scale deployment problems involving only sea 
transportation assets have been made by Buvik [Ref. S] and Lima [Ref. 9] using set 
covering models with dynamic column generation. 

Another initiative which is being pursued by USTRANSCOM to solve the deploy- 
ment problem is the development of the Deployment Analysis Prototype (DAP). The 
DAP is a deployment modeling system developed principally as a technology demon- 
stration implementable entirely on a 80286-based personal computer using many of the 
same inputs as the JDS. This system is used as a base for the development of modeling 
techniques which could be used as stand-alone tools for deployment planning as well as 
a testbed for the development of innovative planning techniques for the JDS. Of specific 
interest for this thesis will be the continuous variable model designed to solve the de- 
ployment problem for air transportation assets which was developed by Rosenthal 
[Ref. 10J to succeed the original air transportation allocation model in the DAP. 

The nature of deployment modeling on a large scale poses many challenges when 
optimization methodology is used. This thesis develops prototypic planning tools to aid 
USTRANSCOM with the allocation of air transportation assets in an effective and 
timely manner. Formulations are developed using the same problem format, but with 
a different modeling approach, of those previously laid out by Rosenthal for the allo- 
cation of air assets. The purpose for developing these new models is to enhance system 
solution-time performance, model resolution, flexibility, and compatibility with other 
modeling programs. Three optimization models and an iterative solution algorithm de- 
signed to improve the DAP are provided. The models and algorithms form a collection 
of ideas related to the deployment problem as it relates to air transportation. 

Chapter 2 of this thesis provides a general discussion of the aircraft deployment 
problem as well as the derivation and description of the optimization models and algo- 
rithm. The first model is based on the weighted distribution problem, more commonly 
known as the generalized transportation problem [Ref. 11]. As in earlier models, this 
model uses continuous variables to represent the assignment of movement requirements 
to aircraft and the assignment of aircraft to missions. The convention of using contin- 
uous variables is justified by the continuous nature of air operations, the aggregate na- 
ture of available data and the level of detail required for the preliminary screening of 
potential operation plans. Unlike the earlier models, this model provides information 
on the movement requirements which are not deliverable at the opportune time for the 
operation plan being considered. This information is particularly important for deciding 
whether sufficient deliveries of particular movement requirements will occur and thus 
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allow the military operation to proceed. Additionally, this model lays the foundation 
for developing the other models. The second model, a transformation of the generalized 
transportation formulation, significantly reduces the number of variables when move- 
ment requirements exist across most route, opportune delivery period, and cargo type 
combinations. This model exhibits a more stable model generation and solution struc- 
ture than the generalized transportation formulation when solved as a linear program. 

A solution algorithm which searches for and then solves subproblems over inde- 
pendent time windows of the original time horizon is described. In many instances, this 
technique may ofTer the ability to model extended time horizons and solve larger in- 
stances of problems. 

The DAP, in its current form, uses a post-solution processor which requires integral 
aircraft solutions. In order to reform the linear programming solution to integers as 
required by this post-solution processor, the resulting nonintegral solutions of the con- 
tinuous variable models are converted to integer solutions. This third model uses a 
greedy rounding model to allocate fractional values making use of any remaining aircraft 
supplies. 

Test results of typical problems are provided in Chapter 3. A typical problem size 
involves nine routes. 80 movement requirements distributed across two cargo classes 
involving 200,000 short tons of freight. 250 aircraft per time period using four different 
aircraft types, with 12 time periods in the time horizon. A summary of the thesis is 
provided in Chapter 4. All of the models and techniques described in this thesis are 
programmed in the General Algebraic Modeling System (GAMS) [Ref. 12] using a S02S6 
based personal computer. GAMS was used because it facilitates rapid changes in the 
model constraints to reflect unique situations for a particular military operation. 
GAMS programs for the models described in Chapter 2 are contained in the appendices. 
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II. PROBLEM FORMULATION 



A. GENERAL PROBLEM DESCRIPTION 

The military aircraft deployment problem involves the shipment of military units 
and equipment between ports of embarkation and debarkation within prescribed time 
intervals using an available pool of aircraft. These shipments are assumed to be 
unidirectional, going from one geographic area to another. An example of this would 
be the movement of a rapid deployment force from the east coast of the continental 
United States to Europe. This chapter presents the derivation of a collection of models 
and an algorithm which can be used to evaluate the general suitability of preliminary 
military operation plans which depend on the solution of this problem. The remainder 
of this section provides a general description of the terminology and elements of the 
models. 

For the models used in this study, the planning horizon is a collection of discrete 
non-overlapping time intervals called periods. The interval of a period depends on the 
distances to be travelled and the speed of the aircraft. Any cargo shipped in a period 
will arrive at its destination within that same period. Additionally, the intervals for the 
different periods may be of equal or varied duration. One reason for allowing the du- 
ration of time periods to vary is because of the increasing uncertainty of events and re- 
quirements the farther the time horizon is extended into the future. Consequently, the 
amount of detail required in the model diminishes with time which allows increasing the 
time interval as periods progress. 

The deployment problem involves the transfer of military units and equipment. Any 
collective group of materials or individuals with a common port of embarkation and 
debarkation which is requested to be moved within a specific time period is called a 
movement requirement. The period of requested delivery is called an opportune delivery 
period since it reflects the period of delivery which is deemed to be most beneficial for 
the objectives of the operation plan being considered. In addition to the opportune de- 
livery period, an interval of permissable periods is defined which contains those periods 
for which movement requirements are allowed to be delivered early or late from the 
opportune deliver.’ period. An available-to-load period is also used which indicates the 
earliest possible period for which a particular movement requirement can be loaded on 
an aircraft. This period represents the availability of cargo, while the early delivery 
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period mentioned above represents the earliest delivery period allowable which still 
meets the objectives of the operation. 

All items within a particular movement requirement will have a single homogeneous 
classification with a standard physical unit. For clarity, the definition of the term phys- 
ical units is noted as being quite distinct from the term military units. While the term 
military units refers to an organized body of troops, physical units refers to the fixed 
quantity of measure for items being transferred. The classifications of movement re- 
quirements are typically passengers, bulk cargo, oversized cargo, or any other appropri- 
ate grouping, while the physical units of measure are typically tons for equipment and 
passengers for military units. These movement requirement classes will be called cargo 
for simplicity. 

The aircraft which can be deployed during a particular period are not treated indi- 
vidually, but rather as a pool of aircraft available to perform missions. The assumption 
that deployable assets are pooled without regard to the location of aircraft at the start 
of a period is reasonable because of the continuous nature of aircraft movements and 
the assumption that all embarkation ports are in the same geographic region. This 
simplification significantly reduces the number of variables required in formulating the 
model since the location of each aircraft docs not have to be tracked. Any of the aircraft 
assigned to the pool may be used on successive tours between a single port of 
embarkation and debarkation during that same period. These tours are designated by 
the forward routes which connect the port of embarkation with the port of debarkation. 
Aircraft missions are defined to be a route, delivery period, and aircraft type combination 
which are used to transport various movement requirements. For this study, the aircraft 
delivery periods are called mission periods. Since the deployment problem assumes a 
unidirectional flow of military units and equipment, aircraft are considered to return 
empty from ports of debarkation. Aircraft are classified by aircraft-type and have asso- 
ciated movement requirements which are compatible with the aircraft design. The air- 
craft are assigned from the pool based on the relative importance of the aircraft. One 
reason for assigning aircraft usage priorities is that certain aircraft are better suited to 
particular military operations. 

The aircraft deployment models of this study use formulations with a channel flow 
concept similar to that originally proposed in the SCOPE model [Ref. 3]. Channels or 
routes are defined only on possible pairs of embarkation and debarkation ports where 
movement requirements exists. The capacity to transport quantities of a requirement 
over a route using a specific aircraft type is calculable. This route capacity is a function 



6 



of many factors such as the expected aircraft speed, the length of the interval of time in 
the period, the distance between ports of embarkation and debarkation, the loading ca- 
pacity of the aircraft, the expected layover time at ports, and the compatibility of aircraft 
types with particular ports. The route capacity is an input parameter available to the 
model. 

The model formulations employ continuous variables to represent shipments of 
movement requirements across routes. A shipment assigned to a mission is defined to 
be delivered entirely within the period of the mission. The quantity of cargo shipped per 
period on a route using a particular aircraft is determined by its route capacity. As such, 
shipments are based on the rate of transfer specified by the route capacity during a pe- 
riod and not discrete plane-loads of cargo. The amount of cargo transferred is limited 
by the total supply of aircraft available. The use of continuous variables is reasonable 
when nontrivial numbers of aircraft are required on any route during a military opera- 
tion. For such cases, the continuous solutions should be sufficiently accurate to provide 
a good preliminary assessment of the logistic supportability of an operation plan. The 
models are required to be able to not only track the closure of all movement require- 
ments to their debarkation ports, but also those portions of movement requirements 
which are delivered earlier or later than their opportune delivery period. This is espe- 
cially important for determining whether a marginal operation plan is supportable. 
Aircraft supply constraints are imposed on the model by using total route capacity re- 
strictions determined by the total number of a particular aircraft type available during 
a period. Various linear side constraints can be included to adapt the model to specific 
operational situations such as capacity limitations at ports or minimum shipment level 
restrictions for movement requirements. 

The objective function consists of costs which are combinations of penalties assigned 
to the different events associated with transporting movement requirements. The cost 
scheme for a model includes relative penalties for (a) the nondelivery of movement re- 
quirements. (b) the early or late delivers' of cargo from the opportune delivers' period, 
and (c) the assignment of cargo and aircraft types. The additive combination of the 
different classes of penalties creates a cost scheme which prioritizes all possible ship- 
ments which could occur using the sarious movement requirements, aircraft, and deliv- 
ers - periods. The specific cost of a transfer may not be economically meaningful; 
hosvever. the relative cost of the transfer in relation to the cost of other types of transfers 
is important. The dominant penalties typically are for delivers - deviations from the 
opportune delisery period and nondelivery of movement requirements. The different 
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types of penalty classes are hierarchically scaled so the dominance of particular penalty 
classes is preserved when they are added to lower priority penalty classes. While the 
objective function minimizes the costs derived from these penalties, the general objective 
of the model is to deliver all movement requirements in a timely manner while efficiently 
using available aircraft. 

B. GENERALIZED TRANSPORTATION FORMULATION 

The original formulation by Rosenthal [Ref. 10] for the aircraft deployment problem 
used dynamic equations to model the shipment of movement requirements and early and 
late deliveries aggregated by mission delivery period in order to reduce the problem size. 
However, this model was unable to provide information about the origination of move- 
ment requirements which were delivered early or late and did not consider the possibility 
of an infeasible problem. The models of this thesis take these aspects of the deployment 
problem into consideration with emphasis being place on the size and solvability of 
problems. This section presents a model which solves the aircraft deployment problem 
using a weighted distribution formulation [Ref. 1 1]. 

The weighted distribution formulation is a generalized transportation network which 
weights the variables in a problem with non-unit coefficients. For the aircraft deploy- 
ment problem, these weighting coefficients are multipliers that represent the route ca- 
pacities which regulate the transfer rates of movement requirements. Side constraints 
may be added to the problem to represent specific operational situations. Many of these 
side constraints may be formulated so they extend the network of this Generalized 
Transportation Formulation (GTF). The primary index sets for the model are: 
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embarkation ports. 
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debarkation ports. 
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routes. 


a e A 


aircraft types. 


h e H 


cargo types. 


t,T e P 


periods of time horizon. 



For notational purposes, t and t are dummy indices representing any period in P . The 
location of / and r in a data parameter or variable determines the context for the use of 
either of these dummy indices. The first index set involving P in a data parameter or 
variable always relates to a mission period while the second index set involving P always 



relates to the opportune delivery period. Additionally, the set of routes R , represents 
the cartesian product of embarkation and debarkation ports. 

For simplicity, the data for the problem is expressed using a single measurement 
system in short-tons (stons), stons aircraft, etc., which are used in a specific period. 
Measurement systems using volumetric or other physical units are possible. The data 
for the problem is: 

MREO nll movement requirement of cargo type h in stons to be shipped over route r 
with opportune delivery period r . 

ALP nll available-to-load period for a movement requirement . 

dni* dnh number of periods which a movement requirement may be delivered early 
or late from its opportune delivery period. 

ASUP„, number of aircraft of type a available in mission period t . 

RCAP a „ cargo capacity in stons, aircraft for aircraft type a over route r during 
mission period t . 

PCAP„ port throughput capacity for movement requirements involving cargo type h 
in stons for port i during mission period i . 

Q’ n/ , cost per ston for shipping cargo type h with opportune delivery period r 
over route r using aircraft type a during mission period i . 

C,),, cost per ston for not delivering cargo within a movement requirement hav- 

ing cargo type h with opportune delivery period t required over route r . 

The satisfaction of movement requirements is accomplished using variables associ- 
ated with available aircraft. The delivery of the movement requirements can be identi- 
fied as belonging to two different classes of assignment. The first class of assignment 
represents the direct satisfaction of a movement requirement with aircraft having a 
mission period the same as the opportune delivery period, henceforth to be called direct 
period deliveries. The second class represents the satisfaction of movement requirements 
by aircraft with a different mission period than the opportune delivery period (yielding 
early or late deliveries of cargo) henceforth to be called cross-period deliveries. Both of 
these classes are represented in the variable X ornh . If there are inadequate supplies of 
available aircraft to satisfy all movement requirements, variables associated with 
undelivered cargo are used. The variables for the undelivered movement requirements 
will be called infeasibility variables. These variables bear a sufficiently large cost so there 
use occurs only when there is no other means of deliver}’. The variables for the model 
are: 
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stons of cargo type h with opportune delivery period x shipped over route r 
using aircraft type a during mission period t . 

N rlll undelivered stons of a movement requirement having cargo type 

h with opportune delivery period z required over route r . 

In addition to the primary index sets, the following derived index sets are necessary. 
The set of periods P in the time horizon must contain the set of periods P m for which 
aircraft arc available to perform missions, that is P m a P. The set of periods P in the 
time horizon must also contain the set of periods P c for the opportune delivery periods 
of the movement requirements, that is P„ a P. 

The set L of all combinations of routes, opportune delivery periods, and cargo types 
for movement requirements is: 



L = { R x P 0 x //}. 

The set La L of allowed combinations depends on the existence of a movement re- 
quirement involving a particular route, opportune delivery period, and cargo type com- 
bination. The set S of all combinations of aircraft types available in a given mission 
period is: 



S = [A x P m ). 

The set S a S of allowed combinations depends on whether there are any aircraft avail- 
able during a given mission period. 

The set Q of all possible combinations for transporting movement requirements on 
different aircraft types during mission periods is: 

Q = { (A x P m e S) x (R x P 0 x H e Z)}. 

The set Q a Q is the allowed combinations of A . R , P m , and // such that nonzero 
route capacity exists and a compatible aircraft and cargo combination exist. Addi- 
tionally, the mission period in P m must fall within the interval containing the early devi- 
ation d r . h and late deviation d rrh from the opportune delivery period in P 0 and must also 
be later than the available to load period ALP r . h for a particular movement requirement. 
This set ensures the combination of primary indices occurs only when there is a move- 
ment requirement and the aircraft and route capacity exists to deliver the requirement 
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within mission periods which fall inside the interval of permissible periods. The set Q„, 
of all movement requirements using a particular aircraft type a and mission period i is: 

Qa, = {R X Po X // I R X Po X // X {( aj )} e Q }. 

The set Q nh of all aircraft types available during the mission periods P m with a particular 
route r , opportune delivery period t , and cargo type h is: 

Qnh = i A x p m I A x P m x {(r, r, h)} e Q). 

In addition to the conditions mentioned above, L , S , and Q may contain other implicit 
constraints which are specific to a given problem. The model is formulated as a linear 
program: 



Minimize 



^ v 4- 

/ j ^artzh ' x artzh ‘ 




(a/,/, Q 



(r, t, A)€ L 



Subject To 



V X an JRCAP an < ASUP al 

(r, ?,/:)€ Q jr 



(aj)eS (1.1) 



V, X a rt-h + ^rth ~ Qnh 

(a,i)e Q r . h 



(r. t,/i) e Z (1.2) 



linear side constraints involving X arrh 



(1.3) 



0 ^ ^arr.h (l- -1 ) 

0 < X nh (1.5) 

Equation (1.1) ensures that all movement requirements are satisfied by actual ship- 
ments or infeasibility variables. Equation (1.2) ensures the number of aircraft deployed 
for missions does not exceed the available pool of aircraft. Equation (1.3) represents 
linear side constraints such as those associated with port throughput restrictions. An 
example of such constraints for ports of embarkation is: 
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(1.3') 



Y X amh < PC A P u V 

(a, t , h) 

A similar set of constraints could be constructed for the debarkation ports. Equations 
(1.4) and (1.5) arc standard nonnegativity constraints. 

The first two constraint sets, equations (1.1) and (1.2) form the generalized trans- 
portation network. It is noted that port capacity constraints (1.3') could be included as 
generalized network constraints by expanding the transportation network to a 
transhipment network. This is not pursued here since not all types of side constraints 
can be handled in this manner. A simple three period depiction for the bipartite graph 
of the generalized transportation network is shown in Figure 1. Nodes are divided into 
those related to the supply of aircraft of a particular type available during a mission 
period (5) and those related to the movement requirements ( L ) . The arcs, represented 
by the variables, are directed from 5 to L . Direct period deliver}' arcs are shown with 
solid lines while the cross period delivery arcs are shown with dashed lines. 

The cost scheme for the variables establishes a set of hierarchical penalties based 
on various aspects of the deployment problem. The primary tier of the cost scheme in- 
volves the timeliness of delivery for the three different classes of assignment. There is 
no penalty for the direct delivery of a movement requirement. Cross period delivery 
penalties of PPEX n are assessed dependent upon whether a delivery is late or early and 
upon the amount of the deviation from the opportune delivery period. The penalty 
structure for cross period variables is an increasing function of the earliness and lateness 
of movement requirement delivery. This function imposes a precedence pattern on the 
model which implies a delivery of one day early is better than one day late which is better 
than two days early and so forth. The following relation illustrates this priciple: 

0 = PPEX,, < PPEX,_ x , < PPEX, +1 , < PPEX,_ 2 , ... 

This pattern progresses throughout the entire interval of permissible periods for delivery. 

Two other penalty structures subordinate to the first are imposed, which allow fur- 
ther refinements of the objectives for assigning aircraft to movement requirements. 
These two additional penalty structures form the second and third tiers of the penalty 
scheme and are based on the relative priority to be used for assigning cargo and aircraft 
types to missions. The second tier of the penalty scheme assesses penalties of HPEX h 
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Figure 1. Generalized Transportation Net w oi k 
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based on the assignment priorities for cargo classes involved with cross period deliveries. 
This penalty tier is used to prioritise cargo in the event of an early or late deliver}’. A 
tertiary tier of penalties A PEX„ , assesses penalties for using particular types of aircraft 
for any type of feasible deliver}'. This penalty tier is used because certain types of air- 
craft may be considered to be more suitable to the military operation being considered, 
while other aircraft may be more valuable as a reserve component. 

The hierarchical penalty structure gives a precedence relation for assessing the im- 
portance of variables when higher level penalties are equal. It is noted that each of the 
three penalty tiers are defined over index sets which are independent of the index sets 
of other penalty tiers. Therefore, the individual penalties within a penalty tier are spe- 
cifically related to its index set. Additionally, all subordinate penalties have a sufficiently 
low value so that they do not influence the penalty structure of higher tier penalties, that 
is, lower tier penalties are smaller than the smallest difference of penalty magnitude be- 
tween any two penalties associated with a higher tier. This creates a cost scheme which 
allows the additive combination of priorities without entangling the various objectives 
represented by the precedence relations of the different tiers. The costs associated with 
particular variables is derived form combining penalties as follows: 

C* t7h = PPEX,, + II P EX h + APEX a . 

Additionally, the cost for infeasible deliveries of C-;, h is set larger than the maximum of 
the feasible variable delivery costs in order to ensure that infeasibility variables are used 
only if feasible delivery is impossible. 

C. MOVEMENT REQUIREMENT REALLOCATION MODEL 

The Generalized Transportation Formulation may have a large number of variables. 
A large problem for an 80286-based personal computer might involve twelve time peri- 
ods. four aircraft types, two cargo types, and nine routes. There is the possibility that 
over 10.000 variables could be generated which would most likely make the problem 
intractable to solve using a general linear programming package. For this reason, a 
transformation of the Generalized Transportation Formulation is proposed in order to 
reduce the number of variables. 

In the original formulation, variables are required to track each possible movement 
requirement that could be assigned to various missions. An alternative method for 
tracking assignments is to create separate variable types for the quantity of cargo 
amongst all opportune deliver}' periods assigned to each mission and for the quantity 
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of movement requirements which are delivered in periods other than the opportune de- 
livery period regardless of aircraft type. With this alternative method of variable gener- 
ation. the direct assignment of cross period deliveries of aircraft missions to movement 
requirements is avoided. This limits the size of the problem mentioned above to under 
3700 variables. This new model is equivalent to the Generalized Transportation For- 
mulation with following exception. In the Generalized Transportation Formulation, the 
type of aircraft was determined for the cargo of movement requirements which was de- 
livered early and late as well as that which was delivered on time. In the new model, the 
cargo which is delivered early or late will be reallocated to cargo of some other period's 
movement requirement. Therefore, this late or early cargo is not specifically tracked as 
a separate delivery on a particular aircraft type, but rather as a part of a different 
movement requirement. The information lost is the aircraft type which the early or late 
cargo arrived on. This information is usually irrelevant since it does not matter which 
aircraft the late or early cargo arrived on, but rather which cargo arrived early or late 
and how early or late it was. The transformation is implemented using an aggregation 
of shipments across opportune deliver} - periods for movement requirements using a 
particular route and cargo type delivered during a specific mission period as follows: 
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Additionally, aircraft types are aggregated for each mission period involved with a 
movement requirement as follows: 




a 
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The transformation of X arrrh into these two new variable types is derived by separat- 
ing the possible variables into direct period deliveries and cross period deliveries. The 
direct period deliveries can further be separated into a positive term involving all possible 
cargo shipped on a particular mission and a negative term involving all early or late 
cargo shipped on the mission. Algebraically, the derivation of the transformation for a 
fixed (/-. t, h) is: 
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The final right hand side of this formula has three terms. The first term represents the 
aggregate movement requirements which are transported on a mission using aircraft type 
a in period t . henceforth to be called the mission delivery variable. The second term 
represents the amount from each of the movement requirements not having an 
opportune delivery period of t , which are transported in mission period r . This is a 
reallocation of movement requirements from periods other than t by regrouping them 
with a movement requirement which has opportune delivery period r . The third term 
represents the reallocation from the movement requirement with an opportune delivery 
period of t by regrouping them with movement requirements having opportune delivery 
periods other than r . The variables associated with these last two terms will be called 
the reallocation variables. 

When the infeasibility variable .V,./, and the movement requirement MREQ r . h are in- 
cluded with this transformation, a restructured version of equation (1.1) of the Gener- 
alized Transportation Formulation is obtained. A conservation of flow relation exists 
for a specific movement requirement node involving the set of indices (r, t. h) . Simi- 
larly, X'„ th is substituted into the remaining equations to give commensurate transformed 
versions of the original equations. The new formulation is called the Movement Re- 
quirement Reallocation Model (MRR) The primary index sets for the model are: 
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embarkation ports, 
debarkation ports. 
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routes. 



re R 


routes. 


a e A 


aircraft types. 


he H 


cargo types. 


t, r e P 


periods of time horizon. 



The data for this model uses the same physical units convention as the generalized 
transportation model. The given data for the problem is: 

MREQ nh movement requirement of cargo type h in stons to be shipped over route r 



ALP nh 


with opportune delivery period t . 
available-to-load period for a movement requirement . 


drth ' dnk 


number of periods which a movement requirement may be delivered early 
or late from its opportune delivery period. 


ASUP a , 


number of aircraft of type a available in mission period t . 


RCAP ar , 


cargo capacity in stons aircraft for aircraft type a over route r during 
mission period t . 


PCAP it 


port throughput capacity for movement requirements involving cargo type h 
in stons for port / during mission period i . 


q;,„ 


cost per ston for shipping cargo type /? with opportune delivery period t 
over route r using aircraft type a during mission period t . 


Qnk 


cost per ston of reallocating a mission requirement of cargo type h over 
route r from a opportune delivery period oft to one oft . 


Q* 


cost per ton for not delivering cargo within a movement requirement having 
cargo type h with opportune delivery period t required over route r . 



The variables for the problem are: 



V 

arth 


stons of cargo type h shipped over route r using aircraft type a during 
mission period t . 


^ mi 


stons of a mission requirement with cargo type h over route r with 
opportune delivery period t which is reallocated to mission period i . 




undelivered stons of a movement requirement having cargo type 
h with opportune delivery period t required over route r . 



In addition to the primary index sets, the following derived index sets are necessary. 
The sets P m , P 0 . and 5 are defined the same as in the General Transportation Formu- 
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lation. The set L' of the cartesian product of routes, periods, and cargo types for 
movement requirements is: 



L' = { R x P 0 x //}. 

The set L' c : L' contains the allowed combinations of route, period and cargo indices. 
At least one of two conditions must occur for an index combination to be a member of 
L' . Either a movement requirement involving a route, the opportune deliver}' period, 
and the cargo type exists, or there are aircraft types available to deliver movement re- 
quirements during the opportune delivery period. The last condition does not require a 
nonzero movement requirement to exist since the model allows the reallocation of com- 
patible movement requirements from other opportune delivery periods. 

The set V of all possible combinations of indices for transporting movement re- 
quirements on different aircraft types during the possible mission periods is: 

V — {A x R x P m x H). 

The set I'c V restricts this set by allowing only combinations a, r, t e P, h such that a 
nonzero route capacity exists and there are aircraft available which are compatible with 
the cargo type. The set IV of all possible combinations of mission periods with move- 
ment requirements is: 



IV = { R x P m x P 0 x H }. 

The set IV cz IT of allowed combinations contains only indices which have nonzero route 
capacity for at least one of the available aircraft types of a specified movement require- 
ment. The aircraft type must be compatible with the particular cargo type of H . Ad- 
ditionally, the mission period in P m must fall within the interval containing the early 
deviation d rrh and late deviation d,, h from the opportune delivery period in P 0 and must 
be later than the available loading period ALP r . h . This set ensures that the combination 
of primary indices occurs only if a movement requirement exists and there are compat- 
ible aircraft and route capacity available which could deliver it using mission periods 
which fall within the interval of permissible periods. 

The set V r . h of all aircraft types which can be used for a particular movement re- 
quirement is: 



Kh = {A \ A x {(r,t,/0} e V). 



IS 



The set V c , of all routes and cargo types which can be used with a particular aircraft type 
and mission period combination is: 

V a! = {R x H | R x H x {(a,/)} e V). 

The set If y ; , of all mission periods for a particular movement requirement is: 

b* = [P m I {P m X i(r,r,h)} e lb {t} 6 P 0 ). 

The set IV nh of all opportune delivery periods for a particular route, mission period, and 
cargo type combination is: 

W' r , h = {P 0 I (P 0 x {(r, r,h)} e b « 6 P m }. 

In addition to the conditions mentioned above above, L' , S , V , and IV may contain 
other implicit constraints which are specific to a given problem. The variables within the 
problem are restricted by sets of allowed combinations of primary indices. The variables 
X' cr!h are defined over the derived set V . while the variables T rr/ , are defined over the 
derived set II' . The problem is formulated as a linear program: 
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linear side constraints involving X' arIh and Y rnh 



(2.3) 
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(2.5) 
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ri-h 
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( 2 . 6 ) 



The interpretation of the constraints in the Movement Requirement Reallocation 
Model is similar to that of the constraints in the Generalized Transportation Formu- 
lation. Constraint (2.1) ensures all movement requirements are satisfied by mission de- 
liver}’ variables, using reallocation variables to transfer movement requirements into or 
out of the node which is represented by the constraint. Constraint (2.2) ensures the 
aircraft pool is not exhausted. Figure 2 shows a simple depiction of the network relation 
obtainable from these two constraint sets. The mission delivery arcs are shown with 
solid lines while the reallocation arcs are shown with dashed lines. Equation (2.3) re- 
presents the linear configuration constraints while equations (2.4), (2.5), and (2.6) are 
nonnegativity constraints. 

Since the network involved with the transformed model is no longer a bipartite 
graph directed from S to L' , there is the distinct possibility of illegitimate flows to exist. 
An illegitimate (low exists when nonzero reallocations occur on a directed path of real- 
location arcs connecting three or more mission requirement nodes. In the event the 
ending node is the same as the originating node for the path, a cycle is created. A di- 
rected path implies that a node on the path will have both an arc importing flow into it 
and another arc exporting fiow from it. If an illegitimate flow exists in the model, then 
a movement requirement is shipped with a lower total cost by using successive reallo- 
cations across intermediate time periods. This avoids incurring the larger cost of real- 
locating a movement requirement directly from the opportune delivery period to the 
mission period. Since the resulting shipment by either reallocation process is the same, 
the use of the intermediate transfers violates the intended relative cost structure of the 
model. 

To prevent the occurrence of illegitimate (lows in the model, a restriction on the 
costs for the reallocation variables and the infeasibility variables of the objective func- 
tion is required. Letting the subset { k,l,m } c= P 0 represent any three time periods of the 
set of opportune deliver}' periods with both r and h being fixed, the following restriction 
will be imposed on the costs involving a movement requirement reallocation between 
periods k and m: 
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Figure 2. Movement Requirement Reallocation Network 
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Lemma: With regard to the Movement Requirement Reallocation Model, condition (1) 
is sufficient to prevent the simultaneous importation and exportation of flow on reallo- 
cation arcs at a L' node. 

Proof: Let node / be the L' node with both route r and cargo type /; fixed and with the 
opportune delivery period of / . Let this node contain a directed path between nodes k 
and m . Condition (1) states the cost for using the infeasibility arc associated with node 
/ is less than the total cost of the two directed arcs connecting it with nodes k and m . 
Since the infeasibility arc associated with node / is unbounded and the Movement Re- 
allocation Model is a minimization problem, the infeasibility arc will always be used 
prior to the use of any combination of two or more reallocation arcs through the di- 
rected path containing the node / . Therefore, the simultaneous importation and 
exportation of flow on reallocation arcs can not occur in the Movement Requirement 
Reallocation Model. QED. 



While condition (1) is sufficient to prevent illegitimate flows in the transformed 
model, additional conditions are required to ensure the use of variables in the models 
are valid. The additional conditions are: 
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Condition (2) states the cost of using the reallocation variable for node k must be less 
than the cost of its infeasibility variable. This ensures reallocations of movement re- 
quirements occur before a movement requirement is determined to be undeliverable. It 
is noted that when restrictions (1) and (2) are assembled together, a triangle inequality 
exists separated by the cost associated with the infeasibility variable. However, the tri- 
angle inequality is not sufficient by itself to prevent illegitimate flows since a directed 
path could exist between nodes k, I, m while no variable existed linking nodes k and in 
directly. The addition of the infeasibility cost to the restriction prevents the simultane- 
ous importation and exportation of flow for this type of directed path which would 
otherwise satisfy the triangle inequality. Restriction (3) is implied by the definition of the 
opportune delivery period. For the definition to be meaningful, movement requirements 
must be satisfied within missions periods which are the same as the opportune delivery 
period prior to the reallocation of the requirement to a node with a different opportune 
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delivery period. Consequently, the cost of satisfying the requirement in the same mission 
period as the opportune delivery period must be less than the cost of reallocating it to 
a different opportune delivery period and the associated cost of performing the mission 
during that period. 

The cost scheme in the General Transportation Formulation is convertible to an 
equivalent structure for the Movement Requirement Reallocation Model if cost re- 
striction (1) is satisfied for both models. It is assumed restrictions (2) and (3) already 
exist in a form useable to the General Transportation Formulation. Because the mission 
delivery costs in the models are independent of the reallocation costs, the usage penalty 
for aircraft are exclusively assessed on the mission delivery variables since they represent 
the actual shipment of cargo by aircraft. Both the timeliness of delivery penalties and 
the cargo assignment penalties represent costs exclusively associated with the reallo- 
cation variables. Therefore, the transformed cost relations are: 

dnh = APEN a 

C y n , h = PPE. Y„ + UPEX,,. 

When comparing the two models, it is important to note the distinct difference be- 
tween the sets L of the Generalized Transportation formulation and L' of the Movement 
Requirement Reallocation Model. Both are used to restrict the first constraint set of 
their respective models. There will typically be many more constraints in the constraint 
set of equation (2.1) of the Movement Requirement Reallocation Model. This occurs 
because the index set L' does not eliminate the node ( r , x, h) when a movement re- 
quirement does not exist for those indices. This null movement requirement node must 
still exist since movement requirements from other opportune delivery periods may be 
reallocated to it. This presents a tradeoff between the two models involving the re- 
duction of variables versus the reduction of constraints. It is expected there will be a 
cross-over point for the solvability of instances of problems dependent upon the density 
of the number of nonzero movement requirements in relation to the total possible 
number of movement requirements given the size of the primary index sets. The General 
Transportation Formulation should perform better for low density instances of prob- 
lems. The Movement Requirement Reallocation Model will improve its performance in 
comparison with the General Transportation Model. 
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D. INDEPENDENT TIME-WINDOW SEPARATION ALGORITHM 



For many operation plans involving the aircraft deployment problem, the size of an 
instance of the model may be intractable or have an unacceptably long solution time. 
This section of the thesis explores a specialized solution technique which can be used to 
solve a certain subset of these difficult problems. The subset of problems involves those 
cases where either the early deliver or the late delivery of cargo is unacceptable. When 
the model is restricted to allowing only early or only late deliveries of cargo, a simple 
algorithm can be used to find independent subproblems within the entire problem. 

Assuming a problem under consideration has no gross inadequacies in the available 
supply of aircraft with regard to the total amount of cargo to be shipped, it is anticipated 
that independent subgroupings of periods may exist. These subgroupings or time- 
windows form a subproblem which is independent of the remainder of the model if the 
movement requirements within the time-window can be satisfied entirely using the air- 
craft available in that same time-window. It is usually difficult to tell where these in- 
dependent time-windows would occur for a particular instance of the problem without 
solving the entire instance. However, it is possible to iteratively define and solve sub- 
problems of independent time-windows when either earliness or lateness is not allowed. 

Consider an instance of the deployment problem where early deliveries are disal- 
lowed. A time-window is independent of future periods if no late deliveries occur outside 
of this time-window. This happens in general because the cost structure of the deploy- 
ment models enforce the delivery of movement requirements as close to the opportune 
delivery period as possible. Since movement requirements can not be delivered in peri- 
ods earlier than the time-window and all deliveries to time periods later the time-window 
have higher costs, then all deliveries are made within the time-window if there are suffi- 
cient assets available. 

An iterative algorithm can be used to solve the entire problem which makes use of 
this property. The algorithm starts by initially solving a subproblem containing a time- 
window’ with only the first period of the time horizon. A time-window of length n is 
a) = { /, t + 1, t -f 2. . . . , t + n — 1}, where t + k e P and k = 0. 1, 1 . If all 
deliveries are made within co , the subproblem is labeled as independent, the solution for 
the subproblcm is saved, and a) is revised to include only the next period of the time 
horizon. If not, the solution is disregarded and the next period of the time horizon is 
added to the existing period in co The next subproblem is solved based on the revised 
to . This process of solving the subproblem and then updating co continues until the last 
period of the time horizon has been solved. Thus, a forward directed algorithm is ob- 
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tained which uses an ascending ordering of the periods (earliest to latest) to solve the 
entire problem by progressively solving subproblems over u> . A similar procedure is 
possible when late deliveries are disallowed which iteratively solves subproblems based 
on a descending progression (latest to earliest) through the periods of the time-horizon. 
A simple algorithm outlining the general procedure is shown in Figure 3. 

This algorithm is particularly well suited to instances of the deployment problem 
where the demand for movement requirements fluctuates substantially from period to 
period or is unimodal in nature with respect to time. As the size of the problems at- 
tempted increases, the benefit in solving the smaller independent time-windows of those 
problems increases. Two other important benefits of the algorithm are its capability to 
model extended time horizons for certain instances, and its ability to be restarted from 
advanced positions in the time horizon. The worst case instance for this algorithm 
would have deliveries outside of the opportune delivery periods for all periods in the 
time-horizon. In this case, the algorithm will ultimately solve the entire problem, being 
no worse off for using the algorithm than having pointlessly expended effort solving the 
smaller subproblems. If the problem is still intractable after using this algorithm, there 
is still valuable partial information which was obtained from iterating through the vari- 
ous time-windows of the algorithm. Even if no independent time-windows were found, 
bottlenecks in the problem would most likely be identifiable from the partial information 
or an advanced starting position found for continuing the algorithm following time- 
windows which are nearly independent. Another use for the algorithm is as a tool to 
determine bounds on instances of problems where both early and late deliveries are al- 
lowed. The procedure would be to solve the problem twice, once where earliness was 
disallowed and once where lateness was disallowed. Since either condition is a restriction 
of the original problem, either offers a bound on the unrestricted problem. By studying 
both solutions, valuable information for the total solution can be obtained. 

E. INTEGER ROUNDING MODEL 

In order to make the solutions from the previous models compatible with the exist- 
ing post-solution processor in the DAP, an integer rounding method is proposed. From 
the solution of any of the previous models, the amount of each movement requirement 
to be transported on the various aircraft types over a specific route and mission period 
are determined. The number of aircraft required for each mission can be obtained from 
this solution using either of the following equations: 
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ALGORITHM: Independent Time-Window Separator 

INPUT: Military Aircraft Deployment Problem formulated 

using a Generalized Transportation Formulation 
or Movement Requirement Reallocation Model. 

Periods P will be preordered from earliest to latest 
(ascending order). 

OUTPUT: List of aircraft mission assignments to 

movement requirement. 

{ Initialize: 

IF (no lateness allowed) 

reorder P in descending order 
co = time-window of the first period of the ordered P 
WHILE ( co e P) 

{ WHILE ( MREQ(a)) = 0) co = next period 
IF( CO E P) 

{ SOLVE the subproblem over co 

D(a>) = number of deliveries outside of co 
IF ( D(co) > 0) 

{ co = co 4- next period } 

ELSE 

{ co = next period } 



OUTPUT 



Figure 3. Independent Time Window Separation Algorithm 
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Generalized Transportation Formulation 



Movement Requirement Reallocation Model 



where the variables on the right-hand side are at their optimal values. The value will 
often be fractional lor AREQ m . 

The rounding model of this study provides one method to round the continuous 
solution. This model attempts to induce integer solutions by rounding up fractional 
values of AREO c „ provided there are unused aircraft available of the same aircraft type 
as the variable being rounded. This avoids leaving cargo undelivered when there are 
unassigned aircraft which could deliver it. The underlying principle is to sort the frac- 
tional part of the continuous solution by magnitude from highest to lowest for each 
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aircraft type during a mission period. The fractional parts are then rounded up to the 
value of one. maintaining the ordering of the sort, until supplies for an aircraft type 
during a mission period are exhausted. Any fractional part of a value which cannot be 
rounded up because of aircraft supply restrictions is omitted (implicitly truncated to 
zero) from the solution. This principle is programmable using a variety of methods. 
The Integer Rounding Model which follows applies this principle in a linear programming 
environment. The primary index sets of the Integer Rounding Model are: 

A aircraft types. 

R routes. 

P m period of aircraft missions. 

H cargo types. 

The given data for the problem is: 

ASUP a , number of aircraft type a available to perform missions during period t . 

AREQ,,, continuous solution from previous model expressed as the number of air- 
craft used. 

E an fractional part of AREQ or , 

The variables for the problem are: 

indicator for rounding up a fractional aircraft of type a used on route r for 
a mission during period i. 

In addition to the primary sets the following derived sets are necessary. The set S 
of all possible combinations of aircraft and mission periods such that aircraft are avail- 
able is the same as described in the previous models. The set / of the aircraft types, 
routes, and mission periods such that fractional aircraft solutions exists is: 

/ = { A x R x P m : E an 0 } . 

The problem is formulated as a linear program: 
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Maximize > E ar ,Z art 



(i a,r,t)ef 

Subject To 



E-ar: — ASlP at 'y ^ ( AREQ art E art ) 



(a,t)eS (3.1) 




E ar i 6 { 0 . 1 } 



(3.2) 



The objective function maximizes the total of the fractional parts of the continuous 
aircraft solution which are rounded up using the binary indicator variable Z m . Con- 
straint (3.1) ensures the indicator variable Z cr . does not exceed the remaining supply of 
aircraft. The remaining supply of aircraft is shown in the right hand side of the equation 
as the total supply of available aircraft less the integer portion of the continuous aircraft 
solution previously calculated. Constraint (3.2) is the zero-one restriction on Z„„ . 

The Integer Rounding Model, when solved as a continuous linear program replacing 
equation (3.2) with 0< Z crt < 1 , will always generate integer solutions. To understand 
this from a mathematical programming perspective, the nature of the constraints of 
equation (3.1) must be explored. Letting U„ r , be the unused aircraft represented by the 
right hand side of equation (3.1). consider the matrix form of the relaxed integer model: 



where A is the technical coefficient matrix. Each of the constraints generated from 
equation (3.1) forms a set of ones for each row of the A matrix. Since the indices of a 
and i are fixed for each equation, mutually exclusive sets of ones are generated for each 
constraint across the routes. Therefore, there will only be a single one in any column 
of A, the rest of its elements being zero. Each of the rows of the A matrix are part of a 
set of generalized upper bound (GL’B) constraints. Additionally, a simple upper bound 
of one is applied to each variable. The proof of integrality follows directly from the 
unimodularity of the A matrix and the integrality of both the variable bounds and U, 
which implies any extreme point is integer [Ref. 13]. 

At this point, it must be emphasized that this model does not necessarily provide 
an optimal conversion of the original deployment model solutions to integer. It is an 



max { EZ : AZ < U, 0 < Z < 1 } 
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approximation. However, the approximation should be close to the optimal solution 
when nontrivial numbers of aircraft are required for most routes during a military oper- 
ation. Notice this condition is the same as that previously mentioned for the deployment 
models. Consequently, a reasonably exact measure using a rounded solution should be 
obtained when an instance for one of the deployment model satisfies this condition. 
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III. COMPUTATIONAL RESULTS 



This chapter describes the results of testing representative data on the models and 
the algorithm previously described. The data for the tests were obtained from 
USTRAXSCOM. The particular data set is similar to what would be expected to be 
solved on the DAP. The models were tested using this data on a personal computer with 
similar characteristics as those which the DAP is currently implemented on. The com- 
puter used was an IBM Personal System 2 Model 50 with 640k of addressable memory 
and a 10 Mhz S02S7 math coprocessor. 

The data set involves four aircraft types, nine routes, twelve time periods, and two 
cargo types. The nine routes are spread over three embarkation ports on the east coast 
of the United States and three debarkation ports in Europe. The number of aircraft 
available during any period ranges from 186 to 354. The test problem requires the de- 
livery of S3 different movement requirements involving 219.693 short tons of freight. A 
crude approximation for the number of aircraft trips required for this problem yields 
over 4000 round trips between the different ports of embarkation and debarkation. This 
large number makes the use of a continuous variable deployment model appropriate. 
For this particular data set, the maximum amount of lateness for any movement re- 
quirement is three days, while the maximum amount of earliness is only restricted to 
periods within the time-horizon. 

Considering the number of index combinations for routes, time periods, and cargo 
types, there are 216 different possible movement requirements which could occur for this 
problem. The S3 movement requirements which actually exist out of the possible 216 
were considered to be representative of the relative density of the actual movement re- 
quirements to the total possible movement requirements for a problem. To examine the 
performance of the General Transportation Formulation (GTF) versus the Movement 
Requirement Reallocation model (MRR), the original set of movement requirements 
(data set ORIG) was expanded or reduced to yield a higher and lower density test in- 
stance. The higher density instance (data set HIGH) contained 135 movement require- 
ments involving 225.103 short tons of freight. The lower density instance (data set 
LOW) contained 33 movement requirements involving 106,906 short tons of freight. 
All other data elements remained as in the original data set. In order to test the Inde- 
pendent Time-Window Separation Algorithm (ITWSA), a second comparison run was 
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performed with lateness disallowed in the model. The ITWSA was applied to the MRR 
model. A test without the separation algorithm was also run for comparison. 

1 he solution characteristics for the linear programs of the various tests are shown 
in Table 1, and the total GAMS execution times and the linear program solver times are 
shown in Table 2 and Table 3. Test type number one is for test runs where lateness was 
allowed in the instance and test type number two is for test runs where lateness was 
disallowed. The type of solution technique is indicated by the name of the model. 



Table 1. LINEAR PROGRAM SOLUTION CHARACTERISTICS 



MODEL DATA SE1 


LOW 


ORIG 


HIGH 


GTF-1 








equations 


144 


196 




variables 


742 


1613 




nonzero coefficients 


3607 


7812 


.... 


MRR-1 








equations 


337 


337 


337 


variables 


760 


1077 


1463 


nonzero coefficients 


3325 


4226 


5332 


MRR- 2 








equations 


337 


337 


337 


variables 


665 


837 


1078 


nonzero coefficients 


3040 


3506 


4177 


ITWSA-2 








equations 


347 


347 


347 


variables 


596 


662 


714 


nonzero coefficients 


2807 


2947 


3791 


max eqns. iteration 


57 


57 


57 


max vars. iteration 


103 


119 


119 


max nzcoef. iteration 


481 


521 


521 


NOTATION: 








GTF-1 Generalized Transportation Formulation 


- run 1 


MRR-1 Movement Requirement Reallocation Model 


- run 1 


MRR- 2 Movement Requirement Reallocation Model 


- run 2 


ITWSA-2 Indep. Time-Window Separation Algorithm 


- run 2 



Table 2. TOTAL EXECUTION TIME (minutes) 



MODEL DATA SE1 


LOW 


ORIG 


HIGH 


GTF-1 


9.92 


25.12 


.... 


MRR-1 


10.62 


13.32 


16.42 


MRR-2 


10.18 


11.67 


13.77 


ITWSA-2 


13.03 


13.4S 


14.37 
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Table 3. LINEAR PROGRAM SOLVE TIME (minutes) 



MODEL DATA SE I 


LOW 


ORIG 


high 


GTF-1 


5.00 


17.12 


— 


MRR-1 


5.75 


7.23 


9.45 


MRR- 2 


5.13 


6.22 


7.50 


ITWSA-2 


1.57 


1.73 


1.S7 



The GTF required a consistently lower number of equations in comparison to that 
of the MRR model. This allowed the GTF to perform well for the lower density in- 
stance: however, as the density increased, its size and solution time increased sharply. 
The increase was so severe, that the GAMS program generator was unable to create the 
linear program for the higher density instance because of inadequate memory capacity 
in the computer. This sharp increase in problem size in relation to the number of 
movement requirements demonstrates that this model's ability to solve the deployment 
problem is highly dependent on the density of the movement requirements. This prop- 
erty is not desirable since a variety of instances of the deployment problem may need to 
be solved. 

The MRR model had a slightly larger size for the lower density instance than that 
of the GI F. The growth in model size was very moderate from the lower density in- 
stance to the original and higher density instances which is in sharp contrast to the rapid 
growth found in the GTF. Despite the large disparity in the number of movement re- 
quirements used in these different instances, the MRR model remained stable in both 
size and solution time. This stability is directly related to the constant number of 
equations found in all instances of the problem, as well as the smaller quantity of vari- 
ables and nonzero coefficients required when compared to the GTF. 

The 1TWSA performed quite well for this particular data set. The algorithm was 
able to identify eleven independent time windows in the twelve period problem. The 
accumulated linear program solution times were extremely low. The total execution time 
was stable: however, it was also larger than the nonseparated MRR model solution time. 
While the solution time was small, the total execution time was larger because of the 
additional overhead of generating all of the subproblems. This data set demonstrates the 
usefulness of the algorithm: however, it is not appropriate to think all data sets would 
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perform nearly as well. The algorithm's performance is based on the ability to find in- 
dependent time-windows within a problem. l or problems with few or no independent 
time-windows, the solution time would be expected to increase substantially. Although, 
the maximum problem size would never be any larger than that which would be required 
to solve a problem without the separation algorithm. 

The Integer Rounding Model was used to transform the continuous solutions of 
previous model test runs to integer values. The rounding model performed well giving 
deployment plans which were consistent with the objectives indicated in the the data set. 
The quality results are directly related to the large number of aircraft required across the 
different routes and time periods. Additionally, it is noted that the time required to run 
the Integer Rounding Model was modest. 

The utility or supportability of the operation plan from the perspective of air trans- 
portation is indicated by the number of aircraft required and the amount of early, late, 
and undelivered cargo. For this particular data set. 14,247 stons of cargo were delivered 
one period early and 2.2S4 stons of cargo were delivered one period late. All of the re- 
maining 203.162 stons of cargo were delivered during their opportune deliver}' periods. 
If the periods are considered to be of one day duration, then the shipment of all move- 
ment requirements used 1455 aircraft-days for the continuous variable solution and 1477 
aircraft-days for the rounded solution. The overall utility of the operation plan would 
depend on these numbers being acceptable for the military actions being considered. 

For comparison, approximations for this deployment problem were calculated using 
the ITWSA. When lateness was disallowed, 14247 stons of cargo were delivered one 
period early while 2284 stons went undelivered. When earliness was disallowed. 11195 
stons of cargo were delivered one period late and 5336 stons were delivered two periods 
late. There was no undelivered cargo when earliness was disallowed. Bounds on the 
amount of early, late, and undelivered cargo are easily determined from these solutions. 
Comparisons of the number of aircraft-days required is not as straight-forward. As the 
amount of undelivered cargo increases the number of aircraft-days required decreases, 
but, as the amount of earliness or lateness increases, the number of aircraft-days in- 
creases or decreases depending of the size of the aircraft and the priority of usage for the 
different aircraft types. For the current penalty scheme, the number of aircraft-days re- 
quired decreased for both the no-earliness and no-lateness approximations. 
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IV. SUMMARY 



This thesis has presented techniques which can be applied to the aircraft deployment 
problem. A formulation based on the generalized transportation problem laid the 
foundation for exploring other methods. This formulation, in itself, was shown to be 
potentially unstable in size and solution time within the limitations of the solver and 
computer when solved as a linear program. This instability depended directly on the 
number of movement requirements in a given instance of the problem. Despite this, the 
model may still be useful with other optimization techniques such as decomposition, 
dynamic column generation, or network algorithms. The GAMS modelling program for 
this Generalized Transportation Formulation is given in Appendix A. 

A transformation of the Generalized Transportation Formulation was performed 
by separating the original transportation variables into those which were used for satis- 
fying movement requirements in the time they were required and those which were used 
to satisfy requirements which were early or late. These variables were aggregated to 
produce an efficient reformulation of the former model. This model was found to be 
quite stable and performed admirably on a personal computer. It is the model of choice 
when the problem is solved as a standard linear program. The GAMS modelling pro- 
gram for the Movement Requirement Reallocation model is presented in Appendix B. 

A time window separation algorithm was proposed as an instrument to attack large 
problem instances. The model provides optimal solutions for problems which allow only 
early or only late deliveries of movement requirements. The model can also be used to 
find bounds on the solution for models allowing both earliness and lateness by restricting 
the earliness or lateness variables. The model attempts to find independent groups of 
time periods within a problem so that the entire horizon of time periods does not have 
to be solved in its entirety for any iteration of the algorithm. Rather, a collection of 
subproblems are solved giving a set of solutions which can be assembled into the total 
problem solution. The MRR model in Appendix B is equipped with the Independent 
Time-Window Separation Algorithm. 

Lastly, an integer rounding model was proposed which allocates the fractional part 
of the continuous solution to the residual supply of available aircraft. This model is an 
effective tool to convert continuous variable solutions to integer when there are not 
substantial numbers of low values to be rounded. Although this model has integer 
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variables, it can be solved as a continuous linear program with guaranteed integrality. 
The model provided quality conversions of the continuous solutions of the former 
models for the data sets tested. The GAMS modelling program for the submodule of 
the Integer Rounding Model is presented in Appendix C. 
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APPENDIX A. GAMS PROGRAM FOR THE GENERAL 
TRANSFORMATION FORMULATION 



$ TITLE GENERAL TRANSFORMATION FORMULATION VERSION (89. 2. 17) 
$0NTEXT 

MICHAEL C. PUNTENNEY 

ACKNOWLEDGEMENT: THIS MODEL IS BASED ON THE ORIGINAL AIR 

LIFT ALLOCATION PROBLEM DESCRIBED IN 
THE TRANSCOM LIFT OPTIMIZER VERSION 
(88.6.14) BY RICHARD E. ROSENTHAL. 

$0FFTEXT 

$OFFUPPER OFFSYMXREF OFFSYMLIST 



DATA 



SETS 

Pm all periods in time horizon / / 
A aircraft types / / 
H cargo types / / 
I ports / / 
R routes / /; 
ALIAS(Pm,Po); 



SETS 

0(R,I) mapping of route to embarkation port 

/ / 

D(R,I) mapping of route to debarkation port 

/ / 

C0M(A,H) compatible aircraft cargo pairings 

/ /; 



TABLE 

TABLE 

TABLE 

TABLE 



ASUP(Pm,A) 

MREQ(R,Pm,H) 

ALP(R,Pm,H) 

RCAP(A,R,Pm) 



quantity of available aircraft ; 
movement requirement ; 

ordinate of available to load period for Pm ; 
freight capacity of aircraft types on routes ; 



PARAMETER PCAP(I,Pm) freight throughput limit at ports ; 



SCALARS 

MAXLTE maximum allowable number of late delivery periods; 

MAXERL maximum allowable number of early delivery periods; 
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DERIVED DATA 



Revise RCAP and ASUP if zero aircraft of route capacity are available 
RCAP(A,R,Pm)$(ASUP(Pm,A) EQ 0) = 0 ; 

ASUP( Pm , A) $( SUM(R , RCAP(A,R,Pm)) EQ 0) = 0 ; 



SET Q(A,R,Pm,Po,H) derived index set for allowed X variables; 

Q( A,R,Pm, Po,H) $C0M( A,H) = YES$(( ORD(Pm) - ORD(Po) LE MAXLTE ) 

AND ( ORD(Po) - ORD(Pm) LE MAXERL ) 
AND ( ORD(Pm) GE ALP(R,Po,H) ) 

AND ( MREQ(R,Po,H) GT 0 ) 

AND ( RCAP( A,R,Pm) GT 0 ) ); 



PENALTY STRUCTURE 



PARAMETERS 

APEN(A) 

HPEN(H) 

TPEN(Pm,Po) 

NPEN 



aircraft priorities 
cargo priorities 

time of delivery penalty from Po to Pm 
aircraft infeasibility override penalty ; 



SCALAR 

ASCALE 

HSCALE 

SCALE 



scaling parameter for APEN 
scaling parameter for HPEN 

scaling factor for solution stability /100. 00/; 



Default HPEN(H) assumes cargos entered in ascending priority: 

HPEN(H) = ORD(H) ; 

Note: APEN(A) values specify order in which assets should be used, 
i. e. , highest value signifies most important to conserve. 

Default APEN(A) assumes assets entered in descending priority. 
APEN(A) = CARD(A) - ORD(A) + 1 ; 

Default TPEN(Pm,Po) 

TPEN(Pm.Po) $ ( ORD(Pm) GT ORD(Po) ) = 

2 + (ORD(Pm) - ORD(Po)) / CARD (Pm) ; 

TPEN(Pm,Po) $ ( ORD(Pm) LT ORD(Po) ) = 

2 + (ORD(Po) - ORD(Pm)) / CARD (Pm) -1/(2* CARD(Pm)); 
TPEN(Pm.Po) $ ( ORD(Pm) EQ ORD(Po) ) = 0 ; 

Determine penalty scale magnitudes ; 

HSCALE = ( 0.9 * l/( 2*CARD( Pm) ) ) / CARD(H) ; 

ASCALE =0.9 * HSCALE / CARD(A) ; 



Scale APEN and HPEN to fit penalty hierarchy 



APEN(A) = ASCALE * APEN(A) ; 

HPEN(H) = HSCALE * HPEN(H) ; 

* Default NPEN greater than any penalty associated with Y and less 

* than twice the minimum penalty associated with Y. 

NPEN = 3.5 ; 



* PROBLEM STRUCTURE 



VARIABLES 

X(A,R,Pm,Po,H) movement requirement assigned to aircraft mission 
N(R,Pm,H) infeasible movement requirements 

COST total penalty costs ; 

POSITIVE VARIABLES X,N; 



EQUATIONS 

BALANCE( R , Po ,H) 
ASSETMAX(A,Pm) 
THRUPUT ( I , Pm ) 
OBJDEF 



maintain freight balance constraints 
observe maximum available assets 
observe throughput limits at ports 
define objective function ; 



OBJDEF . . 

SUM( ( A , R , Pm , Po , H) $Q( A , R , Pm , Po , H) , X(A,R,Pm,Po,H) 

* (APEN(A) + TPEN( Pm , Po) + HPEN(H)$(ORD(Pm) NE ORD(Po)))) 
+ SUM((R,Po,H) , (NPEN + HPEN(H) ) * N(R,Po,H)$MREQ(R,Po,H)) 

=E= COST ; 

BALANCE( R , Po , H) $MREQ( R , Po , H) . . 

SUM((A,Pm)$Q(A,R,Pm,Po,H) , X(A,R,Pm,Po,H) ) 

+ N(R,Po,H) $ MREQ(R,Po,H) 

=G= MREQ(R,Po,H) / SCALE ; 

ASSETMAX(A,Pm) $ ASUP(Pm.A) .. 

SUM( (R,Po,H)$Q(A,R,Pm,Po,H), 

X( A , R , Pm , Po , H) /RCAP( A , R , Pm) ) * SCALE 
=L= ASUP(Pm, A) ; 

THRUPUT( I , Pm) $ ( PCAP(I,Pm) LT +INF ) .. 

SUM(R$0( R , I ) , SUM((A,H,Po)$Q(A,R,Pm,Po,H), X( A ,R , Pm , Po ,H) ) ) 

+ SUM(R$D(R,I), SUM((A,H,Po)$Q(A,R,Pm,Po,H), X(A,R,Pm,Po,H) )) 
=L= PCAP(I,Pm) / SCALE ; 



* 



* MODEL / SOLVER 



MODEL LIFTOPT /ALL/ ; 

OPTION LIMR0W=0 , LIMC0L=0, SOLPRINT=OFF , ITERLIM=5000; 
SOLVE LIFTOPT USING LP MINIMIZING COST ; 



3S 



PARAMETER XA(R,Pm,H,A) display aircraft mission cargo assignment 

XD(R,Pm,Po,H) display late or early movement requirements 
PLANES(R , Pm , A) display aircraft assignment by mission; 

XA(R,Pm,H,A) = SUM(Po, X. L(A,R,Pm,Po,H)) * SCALE ; 

XD(R,Pm,Po,H)$(ORD(Pm) NE ORD(Po)) = 

SUM( A, X. L(A,R,Pm,Po,H)) * SCALE ; 

N. L(R,Pm,H) = N. L(R,Pm,H) * SCALE ; 

PLANES(R,Pm,A) $ RCAP(A,R,Pm) = 

SUM( H$COM( A,H) , XA(R,Pm,H,A) ) / RCAP(A,R,Pm) ; 

OPTION XA: 2: 3: 1; 

OPTION XD: 2: 3: 1; 

OPTION PLANES: 2: 2: 1; 

DISPLAY XA, XD, N. L, PLANES; 
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APPENDIX B. GAMS PROGRAM FOR THE MOVEMENT 
REQUIREMENT REALLOCATION MODEL 

USING THE INDEPENDENT TIME WINDOW SEPARATION ALGORITHM 



$TITLE MOVEMENT REQUIREMENT REALLOCATION TRANF ORM AT I ON (II) 

* WITH INDEPENDENT TIME WINDOW SEPARATION ALGORITHM 

* ( MAINFRAME VERSION - LOOPING SOLVE TECHNIQUE ) 

* VERSION (88. 7. 26), MODIFIED ( 89. 2. 17 ) 

$ONTEXT 

MICHAEL C. PUNTENNEY 

ACKNOWLEDGEMENT: THE MRR TRANSFORMATION IS BASED 

ON THE ORIGINAL AIR LIFT ALLOCATION 
PROBLEM DESCRIBED IN THE TRANSCOM 
LIFT OPTIMIZER VERSION (88.6. 14), 

BY RICHARD E. ROSENTHAL. THE ITWS 
ALGORITHM IS BASED ON VERSION (87.12.23). 

$OFFTEXT 

$OFFUPPER OFFSYMXREF OFFSYMLIST 



DATA 



SETS 

Pm all periods in time horizon / / 

Tm(Pm) window of periods for current solve iteration 

TITER(Pm) iteration counter ; 

ALIAS(Pm,Po,PP,PL); 

ALIAS(Tm,To); 



SETS 

A aircraft types 

H cargo types 

I ports 

R routes 



/ / 

/ / 

/ / 

/ /; 



SETS 

0(R,I) mapping of route to embarkation port 

/ / 

D(R,I) mapping of route to debarkation port 

/ / 

COM(A,H) compatible aircraft cargo pairings 

/ /; 
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TABLE 

TABLE 

TABLE 

TABLE 



ASUP(Pm,A) 

MREQ(R,Pm,H) 

ALP(R,Pm,H) 

RCAP(A,R,Pm) 



quantity of available aircraft ; 
movement requirement ; 

ordinate of available to load period for Pm ; 
freight capacity of aircraft types on routes ; 



PARAMETER PCAP(I,Pm) freight throughput limit at ports ; 



SCALARS 

MAXLTE maximum allowable number of late delivery periods; 

MAXERL maximum allowable number of early delivery periods; 



* DERIVED DATA 

* 



* Revise RCAP and ASUP if zero aircraft of route capacity are available 
RCAP(A,R,Pm)$(ASUP(Pm,A) EQ 0) = 0 ; 

ASUP( Pm, A) $( SUM(R, RCAP(A,R,Pm) ) EQ 0) = 0 ; 



SET W(R,Pm,Po,H) derived index set for allowed Y variables; 

W(R,Pm,Po,H) = YES $(( ORD(Pm) NE ORD(Po) ) 

AND ( ORD(Pm) - ORD(Po) LE MAXLTE ) 

AND ( ORD(Po) - ORD(Pm) LE MAXERL ) 

AND ( ORD(Pm) GE ALP(R,Po,H)) 

AND ( MREQ(R,Po,H) GT 0 ) 

AND ( SUM( A$C0M(A,H) , RCAP(A,R,Pm) ) GT 0 ) ); 



* PENALTY STRUCTURE 



PARAMETERS 

APEN(A) 

HPEN(H) 

TPEN(Pm.Po) 

NPEN 



aircraft priorities 
cargo priorities 

time of delivery penalty from Po to Pm 
aircraft infeasibility override penalty ; 



SCALAR 

ASCALE 

HSCALE 

SCALE 

TOLER 



scaling parameter for APEN 
scaling parameter for HPEN 

scaling factor for solution stability /100.00/ 
zero tolerance for defininig delay / . 01/ ; 



* Default HPEN(H) assumes cargos entered in ascending priority: 

HPEN(H) = ORD(H) ; 

* Note: APEN(A) values specify order in which assets should be used, 

* i. e. , highest value signifies most important to conserve. 

* Default APEN(A) assumes assets entered in descending priority. 

APEN(A) = CARD(A) - ORD(A) + 1 ; 
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* Default TPEN(Pm,Po) 

TPEN(Pm,Po) $ ( ORD(Pm) GT ORD(Po) ) = 

2 + (ORD(Pm) - ORD(Po)) / CARD(Pm) ; 

TPEN( Pm , Po) $ ( ORD(Pm) LT ORD(Po) ) = 

2 + (ORD(Po) - ORD(Pm)) / CARD(Pm) -1/(2* CARD(Pm)); 

* Determine penalty scale magnitudes ; 

HSCALE = ( 0.9 * 1/ ( 2*CARD(Pm) ) ) / CARD(H) ; 

ASCALE = 0. 9 * HSCALE / CARD(A) ; 

* Scale APEN and HPEN to fit penalty hierarchy 

APEN(A) = ASCALE * APEN(A) ; 

HPEN(H) = HSCALE * HPEN(H) ; 

* Default NPEN greater than any penalty associated with Y and less 

* than twice the minimum penalty associated with Y. 

NPEN = 3.5 ; 



tV 

* PROBLEM STRUCTURE 

* 



VARIABLES 

X(A,R,Pm,H) 
Y(R,Pm,Po,H) 
N(R ,Pm,H) 
COST 



cargo assigned by aircraft mission 

late or early delivry of movement requirements 

infeasible movement requirements 

total penalty costs ; 



POSITIVE VARIABLES X, Y, N ; 



EQUATIONS 

BALANCE(R,Pm,H) 
ASSETMAX( A , Pm) 
THRUPUT( I ,Pm) 
OBJDEF 



maintain freight balance constraints 
observe maximum available assets 
observe throughput limits at ports 
define objective function ; 



OBJDEF . . 

SUM( (R ,Tm ,To ,H) , (HPEN(H)+TPEN(Tm,To) ) * Y(R,Tm,To,H)$W(R,Tm,To,H)) 
+ SUM( (R,Tm,H), (NPEN+HPEN(H) ) * N(R,Tm,H)$MREQ(R,Tm,H) ) 

+ SUM((A,R,Tm,H)$COM(A,H) , APEN(A) * X(A,R,Tm,H)$RCAP(A,R,Tm) ) 
=E= COST ; 

BALANCE(R ,Pm,H) $Tm( Pm) .. 

SUM( A$C0M( A ,H) , X(A,R,Pm,H) $ RCAP(A,R,Pm) ) 

- SUM( Po $Tm( Po) , Y(R,Pm,Po,H) $W(R ,Pm,Po,H) ) 

+ SUM( Po $Tm( Po) , Y(R,Po,Pm,H) $W(R,Po,Pm,H) ) 

+ N(R,Pm,H)$MREQ(R,Pm,H) 

=G= MREQ(R,Pm,H) / SCALE ; 

ASSETMAX( A ,Tm) $ASUP(Tm, A) .. 

SUM( (R,H)$C0M(A,H) , X(A,R,Tm,H)$RCAP(A,R,Tm) 

/ RCAP( A,R,Tm) ) * SCALE 

=L= ASUP(Tm, A) ; 
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THRUPUT( I ,Tm)$(PCAP( I ,Tm) LT +INF ) .. 

SUM( R$0(R, I ) , SUM( (A,H)$COM(A,H), X( A,R,Tm,H)$RCAP( A,R,Tm) )) + 

SUM( R$D(R,I), SUM( ( A,H)$COM( A,H) , X(A,R,Tm,H)$RCAP( A,R,Tm) )) 

=L= PCAP( I ,Tm) / SCALE ; 




* MODEL CONTROL 

* 



MODEL LIFTOPT /ALL/ ; 

OPTION LIMR0W=0, LIMC0L=0, SOLPRINT=OFF , ITERLIM=5000; 



PARAMETER XA(R,Pm,H,A) 
YA(R,Pm,Po,H) 
NF(R,Pm,H) 
PLANES(R,Pm,A) 



display 

display 

display 

display 



cargo assigned by aircraft mission 
late of early movement requirements 
infeasible movement requirements 
aircraft assignment by mission ; 



OPTION XA: 2: 3: 1; 
OPTION YA: 2: 3: 1; 
OPTION NF: 2: 2: 1; 
OPTION PLANES: 2: 2: 1; 



* ITWSA CONTROL 



PARAMETER ITWSA determines whether the Independent Time Window ; 

Separation Algorithm is ON or OFF 
( 1=0N; 0=0FF ) 

ITWSA = lj 



SET MREQTEST(PM) terminate zero-MREQ iterations 
BLANKEND(PM) terminate last zero-MREQ solve 
SOLVTEST(PM) terminate solve loop when end test satisfied 
PALL(PM) all Tm if early and late deliveries allowed /T1*T12/; 



SCALARS 

ITER 

INFEAS 

CHECK 

TOTLCOST 

SIGN 

TEND 



iteration counter 

determines infeasibilities for the current iteration 
test for nonzero movement requirement iteration 
total cost of all feasible iterations / 0 / 
backward or forward iteration solve discipline 
ending iteration count number ; 



SIGN = - 1 $ ( MAXLTE EQ 0) + 1$( MAXERL EQ 0) 

- 1$( MAXLTE EQ 0 AND MAXERL EQ 0 ); 
SIGN = SIGN$ITWSA ; 

ITER = CARD( Pm) $( MAXLTE EQ 0) + 1$( MAXERL EQ 0) 

- 1$( MAXLTE EQ 0 AND MAXERL EQ 0 ); 
ITER = ITER$ ITWSA ; 
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TEND = 0 $ ( MAXLTE EQ 0) + (CARD(Pm)+l)$(MAXERL EQ 0) 

- ( CARD (Pm)+1)$( MAXLTE EQ 0 AND MAXERL EQ 0); 
Tm( Pm) $( SIGN EQ -1) = YES $(ORD(Pm) EQ CARD(Pm)); 

0) = PALL(Pm) ; 

1) = YES $(0RD(Pm) 

= YES $(0RD(Pm) 



Tm(Pm)$( SIGN EQ 
Tm(Pm)$(SIGN EQ 
TITER( Pm) $( SIGN EQ -1) 
TITER(Pm) $( SIGN EQ 1) 



= YES $(ORD(Pm) EQ 1); 



EQ 1); 

EQ CARD(Pm) ) ; 



* I TVS A SOLVER LOOP 



SOLVTEST(Pm) = YES; 

LOOP(PP$SOLVTEST(PP) , 

MREQTEST(PM) = YES$(CARD(TM) LE 1) ; 

BLANKEND( TITER) = YES ; 

LOOP( PL$MREQTEST( PL) , 

CHECK = 1$( SUM( (R,Tm,H) , MREQ(R,Tm,H) ) ) ; 

ITER = ITER + SIGN$( CHECK NE 1) ; 

Tm(Pm) = YES$(ORD(Pm) EQ ITER ) ; 

MREQTEST( PM) $( CHECK EQ 1 

OR ITER EQ TEND) = NO ) ; 

BLANKEND( TITER) $( SIGN*( ITER - TEND) GE 0) = NO ; 
LOOP(TITER$BLANKEND( TITER) , 

SOLVE LIFTOPT USING LP MINIMIZING COST ) ; 

MREQTEST(PM) = YES; 

ITER = ITER + SIGN ; 

INFEAS = SUM( (R ,To ,H) , N. L(R,To,H) ) ; 

TOTLCOST = TOTLCOST + COST. L$( TOLER GT INFEAS 

OR ITER EQ TEND) ; 

TITER(Ptn) = YES$(ORD(Pm) EQ ITER ) ; 

Tm(Pm) = Tm( Pm) $( INFEAS GT TOLER) + TITER(Pm) ; 
SOLVTEST( PM) $( SIGN'"( ITER - TEND) GE 0) = NO ) ; 



XA(R,Pm ,H , A) 

YA(R,Pm,Po,H) 

NF(R,Pm,H) 



X. L(A,R,Pm,H) * SCALE ; 

Y. L(R,Pm,Po,H) * SCALE ; 
N. L(R,Pm,H) * SCALE ; 



PLANES(R,Pm,A) $ RCAP(A,R,Pm) = 

SUM( H$COM(A,H) , XA(R,Pm,H,A) ) / RCAP(A,R,Pm) ; 



DISPLAY XA, YA, NF, PLANES, TOTLCOST; 
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APPENDIX C. GAMS PROGRAM FOR THE INTEGER ROUNDING 

MODEL 



$ONTEXT 



INTEGER AIRCRAFT ROUNDING MODEL (88.06.01) 



MICHAEL C. PUNTENNEY - MODIFIED ( 89. 2. 18) 

ACKNOWLEDGEMENT: THIS MODEL IS BASED ON THE ORIGINAL AIR 

LIFT ALLOCATION PROBLEM DESCRIBED IN 
THE TRANSCOM LIFT OPTIMIZER VERSION 
(87.12.23) BY RICHARD E. ROSENTHAL. 

$OFFTEXT 

$OFFUPPER OFFSYMXREF OFFSYMLIST 



* ROUNDING MODULE 



PARAMETER FRAC(A,R,Pm) 
PLAN(R,Pm, A) 



fractional PLANE values 

rounded plane values for deployment plan ; 



FRAC(A,R,Pm) = PLANES(R, Pm, A) - TRUNC ( PLANE S ( R , Pm , A ) ) ; 



VARIABLES 

TOTAL total of fractional values rounded up 

Z(A,R,Pm) rounded fractional values of PLANES ; 



Z. UP( A , R, Pm) = 1; 

Z. LO(A,R,Pm) = 0; 

Z. FX(A,R,Pm)$(PLANES(R,Pm,A) EQ 0) = 0 ; 



EQUATIONS 

OBJF objective function 

ALLOC ATE (Pm, A) allocate fractional values observing supplies ; 



OBJF .. TOTAL =E= SUM( (A,R,Pm) , FRAC(A,R,Pm) * Z(A,R,Pm) 

$(FRAC(A,R,Pm) NE 0)) ; 

ALLOCATE( Pm, A) $ASUP(Pm,A) .. SUM(R, Z(A,R,Pm)) =L= 

ASUP( Pm, A) - SUM(R, PLANES(R,Pm,A) - FRAC(A,R,Pm) ) ; 
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MODEL ROUNDER /OBJF, ALLOCATE/ ; 

OPTION LIMROW=0 , LIMCOL=0, SOLPRINT=OFF; 

SOLVE ROUNDER USING LP MAXIMIZING TOTAL; 



PLAN(R,Pm,A) = TRUNC(PLANES(R,Pra,A) ) + Z. L(A,R,Pm) 

OPTION PLAN: 1: 2: 1 ; 

DISPLAY PLAN ; 
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